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We present an experimental study of hybrid turnstiles with high charging energies in comparison 
to the superconducting gap. The device is modeled with the sequential tunneling approximation. 
The backtunneling effect is shown to limit the amplitude of the gate drive and thereby the maximum 
pumped current of the turnstile. We compare results obtained with sine and square wave drive and 
show how a fast rise time can suppress errors due to leakage current. Quantized current plateaus 
up to 160 pA are demonstrated. 
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Development of a quantum current standard based on 
the discreteness of the electron charge has been one of the 
major goals of metrology for more than twenty years [l| . 
The desired standard would generate the current / = nef 
by transfering n electrons (e) per cycle at frequency /. 
The most precise device of this kind was based on seven 
tunnel junctions in series 0. It reached the relative ac- 
curacy of 10~®, but the maximum current was only on a 
picoampere level. However, closing the quantum metro- 
logical triangle Q, which is one of the prime needs for 
the quantum current standard, requires about 100 pA 
or more. The quest for standards with higher currents 
has involved several candidates, based, e.g., on supercon- 
ducting devices f4|] , surface acoustic waves 5] , and semi- 
conducting quantum dots However, none of them has 
reached metrological accuracy. In this Letter, we study 
the benefits of increased charging energy for the turnstile 
based on a hybrid single-electron transistor (SET), see 
Fig. 1(a), proposed and demonstrated in Ref. [7]. It was 
first realized as a NISIN (N=normal-metal, I=insulator 
and S=superconductor) structure, but the SINIS version 
reported in Ref. [1] promises higher accuracy [§]. 

The operation of the turnstile is based on the super- 
conducting gap A, which expands the stability domains 
of the charge states, see Fig. l(b-c). A small bias voltage 
Vb can be applied since the current is ideally zero when 
\Vh\ < 2A/e. The bias yields a preferred direction of 
tunneling, when the gate voltage Vb{t) is alternated be- 
tween the charge states. The gate voltage is convenient 
to express in normalized units as Vg{t)Cg/e = ng(t) = 
Ugo -I- Agw{t), where Cg, rigo and Ag are the gate capac- 
itance, dc offset and drive amplitude, correspondingly. 
The gate waveform, e.g., sine or square, is expressed with 
w{t). By extending the gate sweep over several charge 
states, current plateaus with n > 1 can be obtained. 
However, the first plateau is optimal for metrology. 

The errors of the turnstile were analyzed theoretically 
in Ref. Q in which, e.g., the tradeoff between the pump- 
ing speed and the error rate was considered. Experimen- 
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FIG. 1: (a) Scanning electron micrograph of a typical turn- 
stile, and a sketch of its measurement circuitry. The bright 
areas are AuPd (normal metal), and the darker parts are Al 
(superconductor), (b) Measured (grey) and simulated (black) 
dc IV curves of sample B. The dc gate voltage was swept dur- 
ing the measurement and hence the envelopes correspond to 
the simulated IV curves with Tig — 0, 0.5. (c) Measured dy- 
namic conductance dl/dVh of sample A as a function of Vb 
and rig. The white areas are stability regions of charge states, 
where the conductance is negligible. They are limited by the 
thresholds for tunneling in the direction preferred by positive 
bias (solid lines) and negative bias (dashed lines) . The pump- 
ing gate signal, when applied, alternates along the thick black 
line. 



tally, the most critical error source appears to be the 
sub-gap leakage, i.e^ residual conductance in the super- 
conducting gap 0, 0, [13] ■ Often the leakage shows up as 
a constant slope at low bias. Then it can be described 
as the ratio ?7„ — i?N / Rn between the normal-state resis- 
tance and the leakage resistance at plateau n. The index 
n = refers to the dc zero-bias leakage. Higher-order 
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tunneling events due to Andreev reflection or Cooper- 
pair-electron cotunneling are known to be leakage pro- 
cesses that affect the operation of the turnstile Q . Possi- 
ble other mechanisms include Cooper-pair breaking [ll[, 
non-equilibrium quasiparticles and nonidealities of the 
tunnel barrier. In Ref. [9] it was shown that the higher- 
order processes are suppressed by the Coulomb blockade. 
Specifically, Andreev reflection vanishes in the regime 
Ug = 1/2 ± Ag^Ar, where ^g,Ar = 1/2 - A/4Ec is the 
threshold amplitude for Andreev reflection at the flrst 
plateau. The charging energy = e'^/2C^ depends on 
the total capacitance of the island, Ce. Note that here 
and in the amplitude threshold discussion below we as- 
sume the bias voltage eVb — A, which minimizes the 
thermal errors f7l|. Together with the threshold for the 
wanted forward tunneling process, Agjt = A/4i?c, ^g.Ar 
yields the requirement Ec > A. In this case, the turnstile 
could be operated at error rates below 10^^, which is the 
prime motivation for studying turnstiles with high charg- 
ing energy. This theoretical result was supported by the 
experiment of Ref. [1], where the leakage rjo was below 
10~^. Yet, typically, the leakage of SINIS structures has 
been of ryo ^ 10~^. 

The first hybrid turnstiles were fabricated with stan- 
dard two-layer electron beam lithography process, and 
the charging energy was of the order of i?c/^B ~ A/fce ~ 
2.5 K ill, iS,]. To increase Ec, we used the trilayer germa- 
nium process [l2| . The samples were metallized in an 
electron-gun evaporator. 

In this Letter, we present data from three samples. A, 
B and C, all with relatively high charging energy, but 
with substantial differencies in i?N- AH measurements 
were performed at temperature of about O.IK. The sam- 
ple parameters are summarized in Table HI All samples 
show leakage that is too strong to be explained by the 
higher-order processes, and that is strongest in the gate- 
open state, unlike expected for Andreev reflection [3] . In 
addition, the leakage parameter rjQ does not seem to de- 
pend on i?N- Albeit counterproductive in general, the 
observable leakage helps us here to demonstrate how the 
leakage at the plateau, rji, can be suppressed from 770 by 
using square- wave drive. 

TABLE I: The parameters of the samples A, B and C. 

An (k») gc/fcB (K) A (/icV) E,/A rjo {xlO~='J 
A 1540 10 250 3.4 8 

B 315 7 210 2.9 4.7 

C 64 8.5 220 3.3 7 



The parameters of the samples were estimated by fit- 
ting their dc IV curves to simulations based on the se- 
quential tunneling model including the cooling effect of 
the SINIS structure [l^. The measured and the simu- 
lated IV curves of sample B are presented in Fig. 1(b). 

The large resistance variation makes this set of sam- 
ples ideal for studying the speed of the turnstile. A high 



pumping frequency can cause two kinds of errors: missed 
tunneling and backtunneling, i.e., tunneling in the di- 
rection against the bias. When the gate signal is ex- 
tended over the backtunneling thresholds, the dashed 
lines in Fig. 1(c), tunneling in the forward direction is 
still preferred, but backtunneling can be statistically sig- 
nificant. With sinusoidal drive and at low frequencies, 
however, tunneling occurs typically before the backtun- 
neling threshold is exceeded. Note that backtunneling 
prevents precise pumping of charge in the case of the 
normal metal SET 

The backtunneling effect can be seen as backbending 
of the current plateaus in Fig. 2(a). The effect is most 
pronounced at low Vb and high Ag. The experimental 
backtunneling slopes are reproduced precisely by our sim- 
ulations based on the sequential tunneling model with 
the sample parameters estimated from independent dc 
measurements. Such an excellent agreement between the 
experimental data and simulations was also observed for 
the other two samples. 




FIG. 2: (a) The current of sample A vs. gate amplitude with 
32 MHz sine wave and at eVb ~ 0.8A. The lines show the 
corresponding simulated currents, (b) The measured plateaus 
of sample C with 300 MHz (solid circles) and 500 MHz (open 
triangles) sinusoidal drives and at Tigo — 0, 0.25 and 0.5. The 
bias voltage was eVh = 1.36 A. 

In Fig. 2(b), we demonstrate with the low-ohmic sam- 
ple C a plateau at 160 pA, which is an order of magni- 
tude higher than in Refs. 0,01 • The experiments cannot 
be reproduced by our constant-temperature simulations, 
since there is already significant heating at the highest 
plateaus. More ideal behaviour is expected at the same 
current level at the first plateau. However, we were lim- 
ited by the maximum frequency of our rf generator, which 
is 500 MHz. 

In Figs. 3(a-b) we show the effect of backtunneling on 
the pumped current in sample B at the first plateau, 
Ugo — 0.5, and eVb = A, with both sine and square wave 
drive. In all cases, there is a flat region between the 
thresholds for forward tunneling (Ag^ft) and backtunnel- 
ing (Ag_bt = 3A/4i?c), shown by the dashed vertical lines. 
Above ^g,bt, backtunneling deteriorates the pumping ac- 
curacy. Since the rise time of the square wave was con- 
stant, about 3 ns, the backtunneling rate did not depend 
on frequency. At 120 MHz, sine and square waves are 
almost identical due to the finite rise time of the square 
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High gate amplitude increases the rate of forward tun- 
neling and thus decreases the effect of missed tunneling. 
How close to Ag.bt the gate amplitude can be set, de- 
pends on temperature. In general, the maximum am- 
plitude is inversely proportional to Ec- In the experi- 
ments below, we used the amplitude half way between 
the thresholds, Ag^^^an = A/2i?c, which is, e.g., about 
0.17 for sample B. The pumping regime is thus drastically 
different from that of the first turnstiles with ^ A and 
i 0.5. 
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FIG. 3: The first plateau of sample B at sine wave (a) and 
square wave (b). The frequencies are 120 MHz (solid trian- 
gles), 64 MHz (open squares), and 32 MHz (solid circles). 
Ag^ft and Ag,bt are marked with the dashed vertical lines, (c) 
The first plateau of sample B as a function of bias voltage 
at 64 MHz at sine wave (dashed line) and square wave (solid 
line), (d) The plateaus of sample C at 16 MHz sine (dashed 
line) and square (solid line) cases. For comparison, we show 
also the sub-gap dc IV curve with gate open (dotted fine). 

At zero temperature, the rate of single-electron tun- 
neling process with electrostatic energy change E^ is 

Fp ~ \J ~ A^/-^N,ie^, where i?N,j is the resistance 
of the junction i. The tunneling rates depend on A, and 
since Ep is roughly proportional to 1/Cs, they depend 
also on the R^^iC-s: time constant. However, if the gate 
amplitude is limited by the backtunneling effect, the for- 
ward tunneling rate at the extreme gate value can be 
expressed as Tp = 7aA/i?N,ie^. The tunnehng rate is 
thus independent of the capacitance. The prefactor 7a 
depends on the exact form of the amplitude limit, e.g., 
7a = a/5/2 for Ag^moan and 7a = V3 for Ag,bt- In Ref. [9[, 
the estimate for the maximum amplitude of the turnstile 
was based on Ag^An which increases as a function of Ec- 
Above Ec — 2A, however, backtunneling sets a more re- 
strictive limit. Hence the maximum current saturates to 
somewhat above 10 pA estimated for E^ = 3A with the 



superconducting gap of aluminum and error rate of 10~^. 

There are two obvious strategies to increase the maxi- 
mum current. Parallelization of turnstiles is an attractive 
option because of the simplicity of the operation of a sin- 
gle device. Another possibility is to find an alternative 
superconductor and to fabricate a device with both high 
A and high Ec. 

An additional benefit of high charging energy is that it 
suppresses the leakage when the gate is not open. In 
Figs. 3(c-d) we demonstrate, as suggested in Ref. [§], 
how a square wave signal can be used to pass the gate- 
open state quickly and hence to decrease the slope of 
the plateau as a function of bias voltage. In Fig. 3(c), 
the plateau of sample B is clearly improved by using the 
square wave. 

The effect of the square wave is more pronounced in 
the case of a fast sample and a low frequency. Then 
the influence of leakage is high in relative terms, and the 
rise times of the sine and square waves of our generator 
are significantly different. This can be seen in Fig. 3(d) 
with sample C and frequency 16 MHz. The slopes with 
sine and square waves are 771 = 9.8 x 10"^ and 771 = 
3.8 X 10^^, correspondingly. Also, both the sine and the 
square wave plateaus extend to higher bias voltage than 
the dc plateau. 

The results of Fig. 3 are traceable to the national low- 
current standard of MIKES |14| , which was used to cali- 
brate the Keithley 6430 current meter. It was then used 
to calibrate the DL Instruments 1211 current amplifier. 
The turnstile could not be measured with the Keithley 
current meter directly, because of its relatively high in- 
put noise. The gain of the DL Instruments 1211 amplifier 
was stable within few parts in 10^^ for several months. 
For each measurement, we determined the input bias of 
the amplifier by comparing different polarities. The un- 
certainty of our measurements was about 10~^. In the 
case of sample B at 64 MHz, also the turnstile reached 
the accuracy of about 10^^. 

To conclude, we have studied the operation of SINIS 
turnstiles with high charging energy. The backtunneling 
effect was shown to limit the drive amplitude and thus the 
maximum current of the turnstile at the highest charging 
energies. Under certain circumstances, a square wave 
drive decreases the influence of sub-gap leakage. We note, 
however, that a square wave tends to heat the island, 
whereas a sine signal can also cool it p^]. Hence, the 
optimum signal may be of some intermediate form. 
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